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1. Introduction 

SU(N) gauge theories with different number of fermions in various representations can play 
an important role in Beyond-Standard-Model phenomenology. It is important to clearly establish 
if a given model is QCD-like, governed by its perturbative fixed point, or conformal with an in- 
frared fixed point (IRFP) at a finite gauge coupling. Identifying the renormalization group (RG) 
structure is not particularly difficult far away from the bottom of the conformal window. Unfortu- 
nately the most interesting case, distinguishing a confining, possibly walking model from a newly 
emerging IRFP, is much more difficult. Such an IRFP likely occurs at strong coupling, requiring 
non-perturbative techniques, and lattice simulations offer the best method to attack the problem. 
Since lattice gauge models are always confining in the extreme strong coupling, identifying a con- 
fining, chirally broken phase does not guarantee QCD-like behavior in the continuum limit. It is 
necessary to start in the weak coupling where direct comparison with perturbation theory is pos- 
sible, and connect it to intermediate/strong coupling. If an IRFP is identified in the process, one 
can conclude that the theory has a conformal phase. If no IRFP is found, one has to show that the 
model exhibits confinement and chiral symmetry breaking. This has to be done before the coupling 
is pushed into extreme strong coupling where universality is lost. While there are well established 
lattice methods to carry out the above program, the calculations are not easy and the emerging 
systematical and statistical errors are frequently much larger than comparable QCD calculations. 
Just as in QCD, it is essential to use different methods and consider different lattice discretizations 
to control the errors when one attempts to distinguish walking and conformal theories and predict 
their critical behavior. 

In this work I describe the application of the 2-lattice matching Monte Carlo renormalization 
group (MCRG) method to calculate the bare step scaling function and connect the weak and strong 
coupling regions[[T]]. This approach has the advantage to the frequently used Schroedinger func- 
tional (SF) method [§-|8|] that it has a free parameter that can be tuned for optimization and it uses 
several observables that provide consistency checks. I apply the method to SU(3) gauge models 
with various number of fundamental fermions. In the Nf = pure gauge case I verify that MCRG 
gives consistent results compared to other methods. The Nf = 16 flavor system clearly shows the 
existence of an IRFP with scaling dimension that is consistent with the perturbative one. In the 
Nf = 12 model I am unable to draw a definite conclusion. The calculation has not identified an 
IRFP and I have not been able to push far enough into the strong coupling to establish a confin- 
ing phase connected to the weak coupling region either. Nf = 12 is very close to the confining - 
conformal boundary and will require a more careful investigation, possibly with a different action 
or renormalization group transformation. 

2. The 2-lattice matching MCRG method 

The 2-lattice matching MCRG method identifies pairs of couplings (j3,/3') corresponding to a 
factor of s change in the lattice spacing, a(j3) = a(fi')/s. A sequence of these pairs can be used to 
predict the bare step scaling function and consequently the running coupling in theories governed 
by one relevant operator. 
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MCRG methods explore the phases and fixed point structure of lattice systems using a real 
space RG block transformation. The block variables are defined as a local average of the original 
lattice variables and a scale s transformation decreases the lattice size by the factor s. By integrat- 
ing out the original variables while keeping the block variables fixed, one removes the ultraviolet 
fluctuations below the length scale sa. The RG transformation does not change the physical cor- 
relation length, but the lattice correlation length after n blocking steps is = s^"^°\ The RG 
can have fixed points only when ^=00 (critical) or £ = (trivial). Flow lines starting near the 
% = 00 critical surface approach the fixed point in the irrelevant directions but flow away in the 
relevant one. After a few RG steps the irrelevant operators die out and the flow follows the unique 
renormalized trajectory (RT), independent of the original couplings. If the flow lines from a pair 
of couplings (j3,j3') end up at the same point along the RT but one requires one less step to do so, 
the lattice spacings at j8 and j3' must differ by a factor of s. In order to identify a pair of couplings 
(j6, j6') with = B,/s we have to show that after n and (n— 1) blocking steps their blocked actions 
are identical. Fortunately one does not need to calculate the blocked actions, it is sufficient to show 
that the expectation values of every operator measured on configurations generated with one or 
the other action are identical. Furthermore it is possible to create a configuration ensemble with 
Boltzman weight of an RG blocked action by generating an ensemble with the original action and 
blocking the configurations themselves This suggests the following procedure for the 2-lattice 
matching: 

1. Generate a configuration ensemble of size L d with action 5(j8). Block each configuration n 
times and measure a set of expectation values on the resulting (L/s n ) d set. 

2. Generate configurations of size (L/s) d with action S(fi'), where /3' is a trial coupling. Block 
each configuration n — \ times and measure the same expectation values on the resulting 
{L/s n ) d set. Compare the results with that obtained in step 1. and tune the coupling j8' such 
that the expectation values agree. 

There is considerable freedom in defining block variables. In this work I use a scale s = 2 
transformation 

Vn.il = PfOj[(l - a)U njtl U n+t i,n + T L U n,vU n+v ^U n+ ^ +v ^Uj 1+2fl J , (2.1) 

where Proj indicates projection to SU(3). The role of the parameter a is to optimize the block 
transformation. While the critical surface of a system is well defined, the location of the fixed 
point itself is not physical, it can be changed by changing the RG transformation. It is important to 
optimize the blocking so its FP and RT can be reached in a few RG steps. The optimal blocking is 
characterized by 

1. Consistent matching between the different operators: along the RT all expectation values 
should agree on the matched configuration sets. Any deviation is a measure that the RT has 
not been reached. 

2. Consecutive blocking steps should give the same matching condition. 
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Figure 1: Left panel: The matching of the plaquette for pure gauge SU(3) theory. Right panel: The bare 
step scaling function Sb(fl;s = 2) for the pure gauge SU(3) system as predicted by different methods. 



3. Simulation results 

The following numerical results were obtained with plaquette gauge action and nHYP smeared 
[ ]T0[ pT] ] staggered fermions. I have generated 2-300 configurations separated by 10 molecular 



dynamics steps at each coupling value. The blocking transformation is given in Eq. 2.1 and I have 
used 3 operators in the matching [jy]. 

3.1 Nf = pure gauge SU(3) gauge model 



The bare step scaling function of the pure gauge SU(3) model was studied in Refs. Q12| ] with 
the 2-lattice matching method. I have repeated some of those calculations with a different block 
transformation and extended them to larger volumes and statistics. Where they overlap, the results 
I present below are consistent with the original calculations. This section mainly serves as a test of 
the method. 

Figure [I] illustrates the 2-lattice MCRG. The left panel shows the matching of the plaquette 
with the blocking parameter a = 0.65 RG transformation. 32 4 volume simulations were done at 
/3 =7.0, and the magenta, blue and red symbols give the values of the blocked plaquette after 
2, 3 and 4 blocking steps. The solid curves interpolate the plaquette values, measured at many 
couplings on 16 4 volumes, after 1, 2 and 3 blocking steps. The 32 4 data match the 16 4 values at 
/3' = 6.49 for all blocking levels. The final blocked volume is 2 4 , but finite size effects are minimal 
as one always compares observables on the same volume. 

The matching can be repeated with different operators and RG transformations. The differ- 
ence between the matched couplings, A/3 = /3 — /3', becomes less and less dependent on the block- 
ing parameter as the blocking level increases. The step scaling function is defined at nt = °° as 
Sb{P) = lim nb ^ooA/3. At finite ny I approximate s&(j3) as A/3 at the optimal blocking parameter, 
defined as where A/3 does not change during the last two blocking steps. This definition clearly has 
systematical errors that can be reduced or eliminated only if the calculation is repeated on larger 
volumes. 
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In the right panel of Figure [j] I compare the MCRG results for the bare step scaling function 
with predictions from other methods, including the 1-loop perturbative prediction, {s = 2) = 
0.59 (dashed line). Note that I show errors only for the MCRG results. Errors from the 16 4 — ► 8 4 
matching are considerably larger than from 32 4 — > 16 4 ; this is mainly due to systematic effects. As 
for the other determinations, I consider the Sommer parameter r [13, [14]], the critical temperature 
T c [15], as well as the SF calculation [^] and a recently proposed alternative based on Wilson loop 
matching (WL) [16] for the running coupling. 

In the scaling regime the different approaches should give the same prediction for Sb(fi). It is 
very satisfying to see the agreement between MCRG, vq and T c even at relatively strong couplings. 
In the range j8 € (6.0,7.0) the predicted values differ considerably from the 2-loop perturbative 
results, but for j8 > 7.0 both the SF, WL and MCRG methods approach s^ en) . The 2-lattice MCRG 
matching method is competitive with other methods in determining the running coupling of asymp- 
totically free theories. 



3.2 N f =16 flavor SU(3) model 

I have used the procedure outlined in the previous section to study the 16-flavor SU(3) model. 
Perturbation theory predicts the existence of an IRFP in the chiral limit at a fairly weak gauge 
coupling. Numerical simulations, based on the SF and WL ratio approaches, have verified the 
backward (towards weak coupling) flow of the gauge coupling Jl7| , jig ). Calculations to identify 
the IRFP are cleanest in the massless limit, but in practice it is easier to introduce a small mass in the 
simulations. I used m = 0.01 — 0.02 but verified that none of the observables had any measurable 
dependence on the mass and for all practical purposes the simulations can be considered as in the 
chiral limit. I have generated 16 4 and 8 4 lattices at various values of the gauge coupling to map the 
bare step scaling function. 

The 2-lattice matching method picks up the relevant operator of a FP. However at the IRFP 
of the 16 flavor model there is no relevant operator in the massless limit. As «/, — > oo matching 
looses its meaning completely since all the flow lines converge to the fixed point. With finite nj, 
one is likely to pick up the flow in the least irrelevant direction. Perturbation theory predicts that 
the critical index governing the flow of the gauge coupling at the IRFP is very small for Nf = 16. 
The gauge coupling is nearly marginal and one expects the matching to follow this coupling. The 
left panel of Figure || summarizes the results. The step scaling function is negative for j8 < 7.0 
signaling the existence of an IRFP. At j3 = 7.5 and 8.0 5j(j3) > 0, so this block transformation has 
an IRFP somewhere around j3 « 7.0. Note that the results presented here are different though con- 
sistent with Ref [p as I have extended the statistics and improved the analysis since the publication. 
Fitting S},(fi) with a linear function in the range j8 € (5.8,7.5) predicts the slope 0.036(12). This is 
somewhat larger but not inconsistent with the perturbative value of 0.022. 



3.3 Nf = 12 flavor SU(3) model 

The Nf = 12 system is very close to the opening of the conformal phase. Calculations with SF 
using unimproved staggered fermions identified an IRFP, but spectral measurements might favor a 
QCD-like interpretation [Q |, [18]-[2(J. I have used the 2-lattice matching MCRG to connect the 



perturbative weak coupling phase to the strong coupling, looking for an IRFP or the emergence of 
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Figure 2: The bare step scaling function st,(j5;s = 2). Left panel: Nf = 16; right panel: Nf = 12. 



a confining, chirally broken phase. The calculation itself is very similar to the Nf = 16 case. I used 
16 4 and 8 4 lattices and simulated nHYP staggered fermions at m = 0.01. 

The results for the bare step scaling function is shown on the right panel of Figure ||[ The 
dashed line indicates the perturbative ^ ert (s = 2) =0.155 value. The numerical results approach 
this for large gauge couplings and show a slight decrease at smaller values. The last data point is 
at j3 = 5.6. For /3 < 5.6 matching with the block transformation of Eq. 2.1 on 16 4 — > 8 4 lattices 



does not give consistent results for the different operators. This does not mean that matching is not 
possible, only that this RG transformation does not reach the FP or its renormalized trajectory in 
3/2 blocking steps. 

The MCRG results indicate that there is no IRFP in the Nf = 12 system for j8 > 5.6. 16 4 
volumes are deconfined at this coupling, so this information is not sufficient to determine if the 
model is conformal or QCD-like. If the system is conformal, here has to be a bulk phase transition 
separating the conformal phase from the confining strong coupling lattice phase at least in the chiral 
limit. I have studied the coupling range from j8 = 5.6 to /3 = 2.5, monitoring the plaquette, blocked 
plaquettes, and several other small Wilson loops. None of these observables showed sign of a 
phase transition. 16 4 lattices remained deconfined for j8 > 4.2 (the smallest coupling I performed 
16 4 simulations) and 8 3 x 16 lattices remained deconfined for j8 > 2.5. All I can say at this point is 
that if the Nf = 12 flavor system with nHYP staggered fermions is conformal, its bulk transition is 
at j3 < 2.5. If it is QCD-like, its lattice spacing at j8 = 2.5 is small enough to deconfine an 8 3 x 16 
volume. On the other hand the plaquette at j6 = 2.5 is around 0.5 (out of 3): this is such a small 
value that smearing becomes problematic, and lattice artifacts likely overwhelm the system. 

It is possible that a different RG transformation has an IRFP at an accessible coupling, or a 
different action, possibly with less efficient smearing, allows the identification of either an IRFP or 
direct connection to a confining phase. 

4. Summary 

Identifying the renormalization group structure of many-fermion theories is a difficult task in 
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models that are near the confining - conformal boundary. In this work I presented a method based 
on Monte Carlo renormalization group that can be used to effectively calculate the bare step scaling 
function and connect the weak and strongly coupled regions. My results for the Nf = pure gauge 
system agree with other determinations for the step scaling function. For the N/= 16 model they 
show the emergence of an IRFP and give the critical exponent of the gauge coupling consistent 
with the perturbative value. The Nf = 12 case is much less clear. With MCRG I could not push 
into strong enough coupling to find a confining phase, neither did the calculation show an IRFR 
Different block transformations or actions could overcome the problems and should be investigated 
in the future. 
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